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We use Perturbative Continuous Unitary Transformations (PCUT) to study the one dimensional Extended 
Ionic Hubbard Model (EIHM) at half-filling in the band insulator region. The extended ionic Hubbard model, 
in addition to the usual ionic Hubbard model, includes an inter-site nearest-neighbor (n.n.) repulsion, V. We 
consider the ionic potential as unperturbed part of the Hamiltonian, while the hopping and interaction (quartic) 
terms are treated as perturbation. We calculate total energy and ionicity in the ground state. Above the ground 
state, (i) we calculate the single particle excitation spectrum by adding an electron or a hole to the system, (ii) 
the coherence-length and spectrum of electron-hole excitation are obtained. Our calculations reveal that for 
V = 0, there are two triplet bound state modes and three singlet modes, two anti-bound states and one bound 
state, while for finite values of V there are four excitonic bound states corresponding to two singlet and two 
triplet modes. The major role of on-site Coulomb repulsion U is to split singlet and triplet collective excitation 
branches, while V tends to pull the singlet branches below the continuum to make them bound states. 

PACS numbers: 71.35.-y, 71.30.+h, 71.10.Li 



I. INTRODUCTION 



The ionic Hubbard model was introduced to describe the 
Neutral-Ionic transition in the charge-transfer organic mixed- 
stack compounds such as TTF-Chloranil^. The transition oc- 
curs by change in pressure or temperature of these materi- 
als^. Since then, in addition to experimental studies^, some 
theoretical investigations for identifying the nature of possi- 
ble phases and optical properties have been conducted on this 
modelS and its extended versions^^. 

On the other hand, in transition metal oxides BaTiOs, 
KNbOs, and KTaOa it was noted that the valence of the transi- 
tion metal ion is ^ +32^2, rather than +4 which is expected if 
a simple ionic picture holds. Therefore it is necessary to con- 
sider effects of strong correlations in these materials, leading 
to the so called ionic Hubbard model (IHM)ii. Although the 
phase transitions at zero^S^^. and finite^Sr™ temperatures for 
this model have been studied by many researchers, the physi- 
cal properties of band-insulating regime of this model has not 
been studied much^l. 

Besides Hubbard-Peierls Hamiltonian^^, the ionic Hubbard 
model is also a candidate for Halogen-bridged transition- 
metal compounds^ known as MX-chains, where M stands for 
transition-metal and X for Halogen. There are many theoret- 
ical and experimental researches on the optical properties of 
these compounds. In NiX-chains, due to large on site interac- 
tion, the system is Mott insulator while PdX- and PtX-chains 
are band insulators^"—. 
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The one dimensional EIHM Hamiltonian is as follows: 

H = y^(-l)V,,, + ^^(4,c,+i,, + /i.c.) 

+t/^ni^^ni,| + (1) 

i i 

where A is one-particle ionic potential, t/4 is nearest- 
neighbor hopping amplitude, U represents on site repulsive 
interaction, and V/2 is nearest-neighbor inter-site interaction. 
c\ g. and Ci^a are electron creation and annihilation operators at 
site i with spin cr, respectively, o- = c\ ^Ci_ij is the electron 
occupation operator and n,; = 71^ ^- When A is dominant 
energy scale, i.e. /S. t,U,V the ground state of this model 
is a band insulator, with electrons and holes and their possible 
(anti)bound-states constituting the excitation spectrum. When 
U dominates, i.e. J7 ^ t, A, F the ground state is Mott in- 
sulator, with doublon and holon excitation. In this paper, we 
investigate the one dimensional EIHM by PCUT method in 
the band insulator region. We take the first term in Eq. ([T]) as 
unperturbed Hamiltonian, and the rest as perturbation. In the 
following, PCUT method is introduced. Then by using this 
method, ground state energy, ionicity, spin and charge gaps, 
one electron and one hole dispersion, as well as spectrum and 
coherence length of the electron-hole excitation for EIHM are 
calculated. 



II. PCUT METHOD 

PCUT or perturbative flow equation method is a high or- 
der perturbative approach for calculating the multi-particle 
excitations and observables in model Hamiltoniansa^— . In 
this method, first the Hamiltonian is mapped, by flow equa- 
tions approach^, to a perturbative effective Hamiltonian. This 
effective Hamiltonian is obtained using a generator— which 
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conserves the number of quasi-particles that are the excita- 
tions of unperturbed part. Then by using Hnked cluster ex- 
pansion theorem22, physical quantities in the thermodynamic 
limit can be calculated employing finite size clusters. To be 
self-contained, here we concisely review PCUT method. 

Assume that the Hamiltonian can be written as a perturba- 
tive problem: 

H = Q + xT, (2) 

where xT is the part of the Hamiltonian, the effect of which 
we would like to consider perturbatively. In order to apply 
PCUT method on the Hamiltonian dU, Q is required to have 
equi-distance spectrum. The difference between two succes- 
sive levels is chosen as the unit of energy. Also T must be of 
the following form: 

+N 

n=-N 

where is a finite number, and r„ is a ladder operator which 
decreases or increases the eigen-values of Q by n units. 

When these conditions are satisfied, the flow equation and 
the quasi-particle conserving generator respectively read, 

deH{e)^[r^{£),H{£)], (4) 

7?,,,W=sigIl(Q,;-g,)ff,,,(^), (5) 

where Qi is the ith eigen value of Q. These equations can be 
used to transform the Hamiltonian ^ to the following effec- 
tive Hamiltonian, 

oo 

Q + CiTfl)T{J^), (6) 

fe=l \T^.\=k,M(r?t)=0 

where \ r^\ ~ k means that rrt is a vector of order k, T{r^) = 
Tmi ■ ■ -Trrik^ M(rrt) = mi + • • • + rrife, and the summa- 
tion is over the components of which can take values 
-A^, -A^+l • ■ ■+N. The C(nl) are known as fractional coef- 
ficients which can be obtained from a set of coupled nonlinear 
first order differential equations. Uhrig and his group, who 
formulated PCUT method, calculated these coefficients up to 
high orders which are available, electronically25. This part of 
the method is quite generic and holds for all Hamiltonians that 
fulfill the above conditions. The salient feature of the PCUT 
effective Hamiltonian is the constraint M{ril) = 0. This con- 
straint causes the effective Hamiltonian to connect only the 
states of Q corresponding to same eigen-values. Therefore, 
interpreting the excitations of Q as quasi-particles (QPs), it 
is sufficient to diagonalize the effective Hamiltonian only in 
sectors with definite number of quasi-particles, e.g. in 0-, 1-, 
and 2-particle sectors. This provides a picture of the low lying 
excitation spectrum of the system. 

In order to calculate the matrix elements of i/cff for an in- 
finite lattice, first it is necessary to decompose the effective 
Hamiltonian to irreducible n-particle operators that are clus- 
ter additive-^. Only for these operators, the matrix elements 



calculated for finite size are valid in thermodynamic limit, i.e. 
there would be no finite size error Therefore in the next step 
one needs to identify minimum cluster size along with suit- 
able form of boundary conditions to calculate required matrix 
elements2i. Note that despite the generic part of the method, 
the choice of cluster depends on the properties of the specific 
Hamiltonian at hand. 

PCUT method has been extensively used in spin systems in 
order to calculate ground state energy, 1 -particle dispersion, 2- 
particle excitation spectrum, and observable&2ti2. However, 
to our knowledge, this method has not been applied much 
to fermionic systemsit. In this paper the aim is to imple- 
ment this method on EIHM Hamiltonian in order to calculate 
ground state energy, 1 -electron and 1-hole dispersions, as well 
as spectrum and coherence-length of electron-hole excitation. 
We use Pade approximation with polynomials of various or- 
ders to enlarge the convergence radius of the obtained pertur- 
bative expansions. 



III. APPLICATION OF PCUT METHOD TO EIHM 

To implement PCUT method on EIHM, we write this 
Hamiltonian as Eq. (|2]l. Now we identify x = t and, 

i.a 

T = i^(c[^c,+i,^ + /i.c.) + u^n,,tn„i 

i.(7 i 
V \ ^ 

+ ^2^"i"-i+l, (8) 



where u = — , and v = — . It is convenient to choose A as 
t t 

energy unit and set A = 1. Note that the popular convention 
for the unit of energy in Hubbard model literature is to take 
t as unit of energy. Moreover, a different sign convention for 
t is used in the literature. This should be noted when com- 
paring our results with other references. It is obvious that at 
half-filling, Q has an equi-distance spectrum and its energy 
difference between two successive levels is 1. Therefore the 
first condition is satisfied. The exact ground state of Q corre- 
sponds to a configuration where all odd sites are doubly oc- 
cupied and all even sites are left empty. Since there are no 
electrons in even sites, nor holes in odd sites, in the sense of 
having no quasi-particles, one can interpret the ground state 
as vacuum. The excitations of Q can be obtained by mov- 
ing an electron from odd site to an even site. Or alternatively 
in terms of quasi-particle picture, when an electron is added 
to even sites and simultaneously a hole is added to odd sites. 
This is a 2-particle excitation that includes two types of par- 
ticles. For obtaining the properties of Hcs in this sector we 
must first discuss 1 -particle sectors corresponding to a single 
electron (hole) added to even (odd) sites. These one particle 
sectors themselves include very important physical informa- 
tion like electron and hole dispersions, which can be directly 
measured in (inverse) photo-emission spectroscopies. 
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Since the eigen-states of Q can be used as basis for H, 
it is possible to write T in the form of Eq. ([3]l. The inter- 
action terms in Eq. ^ compose Tq and with some simple 
calculations we can decompose hopping term to r_i, and 
T+i. Therefore second condition is satisfied with = 1, 
which reflects the fact that hopping is associated with a pair of 
particle-holes. These operators are defined as r„ ~ ^ - r„(i), 
where summation is over lattice points and r„(i) only con- 
nects nearest-neighbor sites, i.e. i and i + 1. The r„(i) are 
defined as: 

V 

ro(i) = uni^frii^i + -Him+i, (9) 



4r+i(z) = |0,t)(t,0|-|t,d)Rt|-U,d)(d,;| 

+ |o,rf)(t,;i-U,t>(rf,o| + |t,;)(d,o| 
+ |o,;)(;,o| - |o,d)(;,t I , ieodd,(io) 

where \a, b) in the above notation means |a),; O \b)i+i, and |0), 
I J,), I t)j and \d) stand for a site occupied with no, spin down, 
spin up, and two electron(s), respectively. The above term is 
essentially a right to left hopping, which has been expanded 
in terms of Hubbard operators for practical purposes. We also 
have r_(_i(z G even) = r^j^(i e odd), r_i(z S even) = 
r+i(i G odd), and r_i(i G odd) = T+i{i G even). In writ- 
ing Eq. (fTOb we have used the following sign convention for 
Fermi operators: When Ci,a or ^ act on a state a phase fac- 
tor (— l)"^' " is produced, where Si^a denotes the total (both 
t and l) number of electrons located in the left side of site i. 
Moreover, at any given site i, spin-^ electron always appears 
to the left side of spin-| electron. As long as one is inter- 
ested in one dimensional lattices the choice of suitable clus- 
ters to find the irreducible interactions is simple while it be- 
comes a tricky task in higher dimensions and for lattices with 
complicated structures. It is obvious from Eqs. ^ and ( fTOl i 
that when T±i{i) acts on two nearest-neighbor sites, it always 
transfers the quasi-particles. But Toii) acts on two nearest- 
neighbor sites with no transfer of QPs. When we calculate 
1 -particle hopping terms and 2-particle interactions in ther- 
modynamic limit based on full graph theory^*^, this property 
of ro(i) in EIHM, in analogy with models involving 4-spin 
interactions'*-, generates larger clusters compared to systems 
with solely 2-spin interactions perturbed around the dimer 



A. Ground state energy and ionicity 

Since i^ofr and Q have the same ground state (unless a 
phase transition takes place)^ and since this state is not de- 
generate, the ground state energy of H per site is given by 
£0 = ^0|iJcff|0), where |0) is the ground state of Q or the 
complex vacuum. We have calculated eo up to order 9 in t, U, 
and V. At this order, from Eqs. (|9]l and (fTOl i it turns out that the 
minimum cluster size is 10 with periodic boundary condition 
to obtain the thermodynamic limit value of eo. Ground state 
energy per site is given in appendix. We have plotted ground 
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FIG. 1 : Ground state energy per site versus t for A = 1 and different 
values of u and v. In each parameter set both plain series and the best 
Pade approximation have been depicted. The plain series expansions 
diverge at some special values of t while the Fade approximant are 
converged in the plotted region. 



state energy per site versus t for different values of dimen- 
sionless interaction parameters u and v in Fig. [T] Note that 
in this figure, each graph represents the best Pade approxima- 
tion. It is seen form Fig. [T|that ground state energy increase 
by increasing U and V. This is physically plausible, as repul- 
sive interactions always increase energy of the system, unless 
a phase transition occurs. Note that since in our convention, 
electron-like bands are represented with negative values of t, 
the negative values of u, v in Fig. [T] correspond to repulsive 
interactions U,V > 0. 

From the analytical expression for ground state energy, we 
can use the Hellman-Feynman theorem to calculate the ionic- 
ity per site which is an important quantity for the ionic Hub- 
bard model and is given byii: 

The ionicity per site versus t for different values of u and v 
has been plotted in Fig. |2] In this figure, only the best Fade 
approximation have been depicted. This figure shows that as 
U is increased the ionicity decreases while increment in V 
increases the ionicity. Therefore the on-site Coulomb inter- 
action U tends to decrease the valence. However, the next 
n.n. term V has the opposite effect. This behavior of V is ex- 
pected, as it stabilizes configurations in which n.n. sites have 
opposite charges with respect to the half-filled background; 
i.e. one is doubly occupied, with its n.n. site being empty. 

If we regard Eq. ( |20b for eo as a series expansion in t, the 
coefficient of can be summed as a geometric series, such 
that the ground state energy per site up to order in band 
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FIG. 2: lonicity per site versus t for A = 1 and different values of 
It and V. In each case two graphs corresponding to the best Pade ap- 
proximation and the plain series have been plotted. The plain series 
expansions diverge at some special values of t. The ionicity per site 
decreases (increases) by increasing in U (V). 



insulator phase is, 

^" = ^ - \ 1 + (12) 

Convergence radius for this expression is obviously given by 
A > |3\//2 — U\. This expression can be interpreted as 
follows: In the band insulating regime, A ;» U^V , when 
t = the ground state belongs to subspace of configurations 
in which odd sites are doubly occupied. By turning on the per- 
turbation t, a second order energy gain can be obtained when 
an electron virtually hops to even neighboring site. From 
above expression and Eq. (fTTT i. we can obtain the ionicity per 
site at the same order as follows, 

/=1-^ S- +0(t'). (13) 

^A + -F-C/)2 

These are perturbative expressions for ground state energy and 
ionicity around atomic limit that are valid for t ^ A-t-|y — J7 
in band insulating phase. In the Mott insulating region, similar 
expressions can be obtained for these quantities at y = in 
the limit i <C J7 — A^^. These expressions are derived based 
on the mapping of the ionic Hubbard model to the Heisenberg 
model^ whose exact ground state is known in ID. 

B. Electron and hole dispersions 

To calculate one-particle dispersion, one must add an elec- 
tron to the even sites of the system, while the corresponding 



hole must be added to odd sites. Since i?off only connects 
those states of Q which have same eigen-values, it can merely 
transfer the electron between even sites and the hole between 
odd sites. Therefore it generates hopping terms of the form 

i(j r—2,4:.,. . . ia 

+ to.hY,'^L'^^'^+ trYhl^hi+r,a+h.C.,{U) 
i(7 r— 2,4,... i(7 

where Eq is the ground state energy and ej^ (^^L) creates an 
electron (hole) with spin a at even (odd) site i. The long range 
hopping amplitudes tr,r = 2,4,... for electron and hole are 
identical up to all orders in perturbation theory. However, as 
for the shortest range hopping amplitude, to (i.e. no hopping) 
only the first order of perturbative expansion for the electrons 
and holes are different and are related by to.e ^ M = to,h + 
where fi = J7/2 -I- is the chemical potential at half-filling 
for Eq. ([T]). Note that if we had started from a particle-hole 
symmetric version of the Hamiltonian, by replacing riia- — > 
n-ia — 1/2 in ?7 and V terms of Eq. ([T]i, we would have fi ~ 
leading to to,e = io,/i- 

The hopping coefficients tr in the thermodynamic limit can 
be calculated on a finite cluster based on linked cluster theo- 
rem^. To calculate these coefficients, first we must determine 
the minimum cluster size which can give us the thermody- 
namic limit results. The minimum cluster size, in addition to 
order of perturbation, depends on the particular hopping co- 
efficient which we would like to calculate. For example, to 
obtain to up to order 9, 13 sites and up to order 10, 15 sites 
are needed. For calculation of 1-2 up to order 9, 10, and 11 one 
needs 13, 13, and 15 sites, respectively. For t4, tg, ig, up to 
order 9,11 sites are sufficient to calculate them in thermody- 
namic limit. These numbers can be inferred from Eqs. (|9]l and 
( fTOl ). We have calculated electron dispersion up to order 9. At 
this order there are five hopping coefficient, to, t2, t^, tg, and 
tg which are independent of the spin of the electron and have 
been reported in appendix. 

Electron and hole dispersions of the effective Hamiltonian 
are given by, 

4 

^e{h){k) = ±(to,e(/i) + 2 ^ t2n cos(2nA:)). (15) 

In Fig. [3] we have plotted the one-particle dispersion corre- 
sponding to t = -1.0, U = 1.4, y = 0, as well as t = -1.0, 
U — \A, V ~ 0.3. Note that to obtain reliable results for 
the value of J7 = 1.4 we have used Pade approximation. The 
charge gap is defined as, 

E,, = Ea{N + l) + Eo{N -l)-2Ea{N). (16) 

For C/ = 0, y = 0, we have a simple band insulator with 
gap magnitude A which has been taken to be unit of energy. 
In Fig.[3]the charge gap in terms of the effective Hamiltonian 
becomes, u]e[k = f ) ~ ^h{k = ^). As can be seen for t ~ 
— 1, U = 1.4, V = 0.0, the gap becomes very small k, 0.014, 
which agrees very well with corresponding DMRG valuoii 
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FIG. 3: Electron and hole dispersions around Fermi surface ^ = 
^ + V for a = -1.0, U = 1.4, V = 0) and {t = -1.0, U = 1.4, 
V — 0.3) at A = 1. We have used Pade[7,2] in both cases. Our 
result for charge gap at {t = —1.0, U = 1.4, V" = 0) is 0.014 which 
is in excellent agreement with DMRG result that is 0.012. 
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FIG. 4: Charge gap versus t for different values of u and v in band 
insulator region. In all cases Pade approximation [p, g] correspond- 
ing to quotient of polynomials of order p, q has been used, which is 
indicated in the legend. Spin and charge gaps are equal, because the 
hopping coefficients, i.e. f„, are independent of spin. 



0.012. Note that when comparing results of Ref.ljjjwith ours, 
one has to appropriately re-scale the Hamiltonian parameters. 
Beyond IJ « 1.4 the band insulating phase terminates and 
one may have different behaviors for spin and charge gaps'^. 
Although our perturbative treatment does not extend beyond 
the band insulating state, extension of the convergence radius 
captures a critical value of U at which the charge gap tends to 
zero. Now at [/ = 1.4, we turn on V . As can be seen in Fig. [3] 
turning on V increases the gap. 

Charge gap versus t for different values of u and v have 
been plotted in Fig. |4] For each value it, v the best Pade fits 
have been shown. For u = — 1.4, u = the Pade fit shown 
in the figure approaches the DMRG valusii as t — > —1. For 
other values of w, w reported in this figure, to ensure the quality 
of results, for each value (u, v) we have shown two Pade fits. 
Comparing the variation of charge gap at w = 0, as can be 
seen by moving from u = —0.9 (red dashed line) to u = — 1.4 
(black solid line), the gap magnitude decreases, which means 
Hubbard U work against the band gap parameter A and tends 
to reduce it. Similarly, for a fixed value of u, turning on the 
n.n. repulsion V = vt increases the band gap. 

Similar to ground state energy we find expressions for to.e, 
and toji, around atomic limit up to order P for V = 0. These 
expressions are as follows: 



to,e — to.h + U 



A A + U 



2 8 A(A - U) 
and for t2 at finite V up to order we obtain: 



O(t^), (17) 



t2 



1 



16A + V-U 



0{t'). 



(18) 



From these expressions we can find the electron and hole 
dispersions for V = Q. The dispersions are as u){k) = 
to + 2t2 cos(2fc). This gives rise to inverse square root density 
of states- a characteristic of ID bands, the width of which is 
controlled by t2- 

In this subsection we obtained some informations about 
EIHM at half-filling by adding one electron and one hole to 
the system, separately. This gave rise to the band disper- 
sion relation of the system. Now let us study the particle-hole 
(charge neutral) excitations of the system. 



C. Spectrum and coherence length of electron-hole excitation 

The first excited state of Q at half-filling is when one elec- 
tron from an odd site with energy — ^ is transfered to an even 
site with energy +-i or in quasi-particle language when one 
electron is added to even sites and one hole is added to odd 
sites, simultaneously. This is a 2-particle excitation which in- 
clude two types of particle, i.e. electron and hole. It is possible 
to discuss the spectrum and coherence length of this electron- 
hole excitation by PCUT method. State of the system in this 
sector can be represented as |i, rrie'-j, rrih) where i is the even 
site and j is the odd site at which electron and hole excitations 
are created, respectively. Also m,, and m/i are magnetic quan- 
tum numbers of electron and hole. To discuss the two-particle 
sector, in addition to quadratic (hopping) terms, one needs to 
include quartic electron-hole terms of the following form in 
the effective Hamiltonian, 
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FIG. 5: Plain series result for electron-hole energy versus K for t = 
—0.5, U = 0.5, and = 0. There are three singlet bound/anti- 
bound state branches and two triplet bound state branches. The two 
singlet anti-bound state modes and two triplet bound state modes are 
degenerate at if = ±-|. 
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FIG. 6: Inverse of coherence length versus K for t = —0.5, 
U — 0.5, and V — 0. For 5*1 mode the coherence length remains 
finite in the whole of range of momenta. But for other modes it be- 
comes infinite at some values of K. Coherence length of 5*2 has a 
discontinuity at K* = 0.467r. 



Since total iJcff preserves the total spin of the electron-hole 
pair, we must use the basis \i,j;S, M), where S is total spin 
of electron-hole pair and M is its magnetic quantum number 
S can be either or 1 corresponding to singlet and triplet to- 
tal spin states. When Hcs acts on states with a particle-hole 
pair, it gives rise to a new configuration in the subspace with a 
particle-hole pair. The processes of mapping the initial state to 
final state can be classified as 2-particle, a mixture of 1- and 2- 
particle, or a mixture of 0-, 1-, and 2-particle process. In order 
to discuss this 2-particle sector, first one has to determine the 
irreducible 2-particle interactions obtained by subtraction of 
0- and 1- particle contributions which by linked cluster theory 
can be employed to calculate ti-^ -j-^ j2 for ^ thermodynamics 
system on a finite clustes^LZ^. Again similar to 1 -particle sec- 
tor, the minimum cluster size can be identified from properties 
of Eqs. (|9|l and (fTOl i. We have calculated these coefficients up 
to order 8 in t, U, V. These coefficients turn out to be solely 
dependent on total spin S and not on magnetic quantum num- 
ber M. 

When the interaction terms are identified, we have only a 2- 
particle problem that must be solved. Such a problem can be 
solved by changing the basis to total particle-hole momentum 
and their relative distance \K,d), or to total momentum and 
relative momentum \K,K') of particle-hole pair—. We use 
\K, d) representation, as in this basis it is possible to explicitly 
calculate electron-hole coherence length defined asii: 

L=%4^^, (19) 

where fd is the probability amplitude for finding the particle- 
hole pair at a relative distance d. In our case d can take only 



odd values as electron are always in even sites and hole in odd 
sites. 

Now that we have the explicit form of the _ffoff in 2-particle 
sector, it remains to diagonalize Hd,d'{K) for a given value 
of center of mass momentum K. Here d, d' indicate relative 
distance. For problems where the two particles are of same 
type (e.g. in Cooper pairing channel), this matrix is a semi- 
infinite^-, but for problems with two types of particles, such as 
our case (i.e. in particle-hole channel), it is an infinite matrix. 
However, when dealing with these matrices on the computer, 
one constructs finite N x N dimensional representations of 
the Hamiltonian for which the eigen-value problem must be 
solved. Next by systematically increasing N, the spectrum at 
N —i' oo limit can be obtained. For those values of K where 
the particle-hole pairs form a bound or anti-bound state, the 
convergence is achieved quickly. This numeric observation is 
physically plausible, because for split-off states, the two parti- 
cles are close to each other and hence the finite size effects be- 
come irrelevant. By determining eigen-vectors of the matrix 
we obtain the fd which can be used to calculate the coherence 
length via Eq. ( fT9] l. 

Upper and lower edges of the particle-hole continuum 
(PHC) is obtained by ignoring the 2-particle interactions and 
diagonalizing i/cff in the basis of \K, K') which gives rise to 
free particle-hole energies. Then the maximum (minimum) 
eigen-value in K' for a fixed value of K in Brillouin zone 
determines the upper (lower) edge of the PHC^** "*"*. 

Plain series result of electron-hole excitation spectrum for 
t -0.5, U ^ 0.5, and y = is depicted in Fig.jS] This 
figure shows that there are three singlet and two triplet modes. 
The singlet channel contains two anti-bound states and one 
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FIG. 7: Plain series result of electron-hole energy versus K for 
t = -0.8, U = 0.0, and V = 0.1. There are two singlet (SI 
and S2) and two triplet (Tl and T2) bound state modes. Each sin- 
glet is nearly degenerate with one triplet which show the effect of on 
site interaction U in separation of singlet and triplet state. 



bound state mode while the triplet modes are bound states. 
The two singlet anti-bound state branches and the two triplet 
branches are degenerate at the edges of the Brillouin zone, 
K = Hence, when longer range Coulomb interaction 

(V) is absent, on-site interaction U stabilizes the triplet ex- 
citations with respect to singlet ones. Note that in the limit 
U ^ 0, a simple band picture would imply the singlet and 
triplet branches to have the same excitation energies. There- 
fore such singlet-triplet splitting is a remarkable effect of short 
range correlations which manifests itself at molecular— as 
well as infinite size systemsi^ where the coiTelation effects 
is believed to play important role. 

We have also plotted the inverse of the coherence length in 
Fig. |6] corresponding to parameter values of Fig. |5] Except 
for 5*1 and 52 branches, the coherence length increases with 
decreasing bounding energy. Coherence lengths of 51 and 
52 decrease with decreasing of anti-bounding energy and 52 
mode has a discontinuity at K* = 0.467r. For 51 mode which 
is anti-bound everywhere in the BZ, the coherence length al- 
ways remains finite. But for other branches who enter the 
PHC at some momentum, the coherence length become infi- 
nite at those particular momenta. 

Plain series results for electron-hole energy and inverse of 
coherence length versus K for t ~ —0.8, U = 0.0, and 
V = 0.1 have been depicted in Figs. [T] and [8] respectively. 
Fig. |7] shows that there are two singlet and two triplet modes 
which all of them are bound states. The fact that the all modes 
are bound state compared with previous case, V = 0, in- 
dicates that the long range part of the Coulomb interaction 
amounts to attraction among electron-hole pairs in both sin- 
glet and triplet channels which in turn forms the so called ex- 



citonic states. Each of singlet branches is nearly degenerate 
with the corresponding triplet branch. This again emphasizes 
that the short range Hubbard interaction U is the major mecha- 
nism behind the splitting between singlet and triplet collective 
mode branches. 

In Fig. [8] we display the behavior for coherence length for 
the same parameter values as in Fig. [T] Coherence length 
of 52 and T2, that are nearly degenerate, monotonically in- 
creases by decreasing K until a certain value K* is reached, 
beyond which the coherence length becomes infinite. But for 
51 and Tl there is always a finite coherence length, indicat- 
ing the corresponding p-h bound state exists everywhere in the 
BZ. 
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FIG. 8: Inverse of the coherence length versus K for t = —0.8, 
U — 0.0, and V = 0.1. Inverse of coherence length is identical for 
52 and T2 and increases monotonically with K while for SI and 
Tl there is a peak not exact at the border. 



Now let us consider a case where both U and V are present. 
In Figs. |9] and [TOl we have plotted the plain series results for 
particle-hole energy and inverse coherence length versus K 
for t = -0.8, U = 0.4, and V 0.3. Again the behavior of 
coherence length and binding energy are in accordance with 
each other. As can be see while due to presence of non-zero 
V, all modes are binding, still on-site correlation U leads to 
singlet-triplet splitting as expected. As can be seen in Fig. |9] 
near the BZ edge, triplets are lowest lying excitations. But for 
smaller momentum transfer K, there are regions where the 
singlet mode will become lowest excited state. 

In the appendix we give analytic expressions for the 
electron-hole energy Ep\i{K) for singlet and triplet modes at 
the zone edge, K = ^, up to order 8. These relations are valid 
for < V. 
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FIG. 9: Plain series result of electron-hole energy versus K for t — 
-0.8, U = 0.4, and V = 0.3. There are two singlet, SI and S2, 
and two triplet, Tl and T2, bound state modes. For iv < 0.16 the 
lowest energy mode is singlet while for K > 0.16 it is a triplet state. 
At the border of Brillouin zone, K — ± ^ , the two singlets and two 
triplets are degenerate. 
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FIG. 10: Inverse of coherence length versus K foi t = —0.8, U — 
0.4, and V = 0.3. There are four graphs correspond to two singlets, 
SI and S2, and two triplets, Tl and T2, modes. At the A' = ±f all 
states have same coherence length that is equal to 1. 



IV. SUMMARY DISCUSSION 

We used perturbative continuous unitary transformation 
method to study the ID extended ionic Hubbard model in the 
band insulating regime. EIHM in addition to the usual ionic 
Hubbard model includes a n.n. neighbor Coulomb repulsion, 
V. We treated the ionic potential as unperturbed part and 
studied the effect of hopping and Coulomb interaction terms, 
perturbatively. The convergence range of the perturbative re- 
sults were extended by using a Pade scheme to fit the results 
and extrapolating them to larger values of interaction param- 
eters U, V which enables us to approach the limit where band 
gap vanishes. Effective Hamiltonian was obtained and diago- 
nalized in zero-, one-particle and particle-hole sectors. Zero- 
particle sector gives us the ground state energy which upon 
using the Hellman-Feynman theorem lead to a perturbative 
expression for the ionicity. tonicity in the absence of U, V 
terms at zero temperature has its maximum value of 1, which 
is a characteristic of band insulators, tonicity decreases by in- 
creasing U, while by increasing V, it tends to increase. There- 
fore the on-site correlations U are useful to remedy the va- 
lence discrepancy in transition metal oxides such as BaTiOa, 
KNbOa, and KTaOg. 

In one-quasi-particle sector we calculated the electron and 
hole dispersions which was achieved by adding one electron 
or a hole to the system. Study of the charge gap atV ^ in- 
dicated that the band gap tends to zero for a critical value of U 
which is in agreement with state of the art DMRG results^^^. 

In the particle-hole sector, the total spin of particle-hole pair 
is a good quantum number which helps to classify particle- 



hole excitations in singlet and triplet channels. The particle- 
hole continuum is characterized by the spectrum of free 
particle-hole pairs. Turning on the interaction gives rise to 
split-off (bound or anti-bound) states in the triplet and singlet 
PH channels. There are in general three/two singlet and two 
triplet branches in the reduced Brillouin zone (— 7r/2 < k < 
tt/2). Given that the bound states in each total spin channel 
are degenerate at k = ±7r/2, the two branches can be thought 
of as a single branch in the larger BZ (corresponding to the 
A = Hmit). The so obtained triplet and singlet correspond 
to a pole in the appropriate susceptibility when the problem is 
viewed from itinerant limit. In this sense we can interpret the 
triplet and singlet split-off states as collective modes. We also 
discussed the coherence length of these collective states. In 
general when the binding (anti-binding) energy of the collec- 
tive states which is defined as the energy difference between 
the energy eigen-value corresponding to the collective state 
and the lower (upper) bound of the PHC decreases the coher- 
ence length which is a measure of PH bound- (anti-bound-) 
state increases. At a critical value of center of mass momen- 
tum K where the collective mode branch enters the contin- 
uum of free particle-hole pairs the coherence length becomes 
infinitely large. Within PCUT when the U term is zero, for 
finite V the singlet and triplet collective modes are almost de- 
generate. The major effect of n.n. Coulomb repulsion V is 
to render the singlet branches into binding collective states by 
pushing them below the PHC, while the major role of U is 
to give rise to remarkable singlet-triplet splitting^. Both sin- 
glet and triplet states which will be bound states in presence 
of V can be interpreted as the excitonic states of the band in- 
sulator under study. Excitonic states are well-known for their 
optical signatures. Therefore the PCUT approach seems to be 
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promising for investigation of excitonic effects in the optical 
spectra^MS. Within the PCUT formulation it is also possi- 
ble to discuss the propagation of two-electron and two-hole 
states which touches upon the problem of superconductivity 
in strongly correlated electron systems. 
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Ground state energy and hopping terms up to order 9 in the exponent summation of t, U, and V. 
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